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Abstract. Although the Gauss-Bonnet term is a topological invariant for general 
relativity, it couples naturally to a quintessence scalar field, modifying gravity at solar 
system scales. Wc determine the solar system constraints due to this term by evaluating 
the post-Newtonian metric for a distributional source. We find a mass dependent, 1/r'' 
correction to the Newtonian potential, and also deviations from the Einstein gravity 
prediction for light-bending. We constrain the parameters of the theory using planetary 
orbits, the Cassini spacecraft data, and a laboratory test of Newton's law, always 
finding extremely tight bounds on the energy associated to the Gauss-Bonnet term. 
We discuss the relevance of these constraints to late-time cosmological acceleration. 
Keywords: dark energy theory, gravity, string theory and cosmology 



1. Introduction 

Supernovae measurements [1] indicate that our universe has entered a phase of late- 
time acceleration. One can question the magnitude of the acceleration and its equation 
of state, although given the concordance of different cosmological data, acceleration 
seems a robust observation (although see for criticisms). Commonly, in order to 
explain this phenomenon one postulates the existence of a minute cosmological constant 
A ~ 10^^^ eV^. This fits the data well and is the most economic explanation in terms of 
parameter (s). However such a tiny value is extremely unnatural from a particle physics 
point of view [3]. Given the theoretical problems of a cosmological constant, one hopes 
that the intriguing phenomenon of acceleration is a window to new observable physics. 
This could be in the matter sector, in the form of dark energy [H [5], or in the gravity 
sector, in the form of a large distance modification of Einstein gravity O [71 [H [9]. 
Scalar field driven dark energy, or quintessence |1] is one of the most popular of the 
former possibilities. However these models have important drawbacks, such as the 
fine tuning of the mass of the quintessence field (which has to be smaller than the 
actual Hubble parameter, Hq ~ 10~^^eV), and stability of radiative corrections from 
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the matter sector [TO] (see however [H]). Modified gravity models have the potential 
to avoid these problems, and can give a more profound explanation of the acceleration. 
However, these are far more difficult to obtain since Einstein's theory is experimentally 
well established [I2], and the required modifications happen at very low (classical) 
energy scales which are (supposed to be) theoretically well understood. Furthermore, 
many apparently successful modified gravity models suffer from instabilities or are 
incompatible with gravity experiments. For example the self-accelerating solutions of 
DGP [H] suffer from perturbative ghosts [I3], and f{R) gravity theories [9] can conflict 
with solar system measurements and present instabilities [H]. 

In this paper we will consider observational constraints on a class of gravity theories 
which feature both dark energy and modified gravity. Specifically, we will examine 
solar system and laboratory constraints resulting from the response of gravity to a 
quintessence-like scalar field, which couples to quadratic order curvature terms such 
as the Gauss-Bonnet term. Such couplings arise naturally [I5], and modify gravity at 
local and cosmological scales [IHl [16]. Although the Gauss-Bonnet invariant shares 
many of the properties of the Einstein-Hilbert term, the resulting theory can have 
substantially different features, see for example [1^. It is a promising candidate for 
a consistent explanation of cosmological acceleration, but as we will show, can also 
produce undesirable effects at solar system scales. 

In particular, we will determine constraints from deviations in planetary orbits 
around the sun, the frequency shift of signals from the Cassini probe, and table- 
top experiments. In contrast to some previous efforts in the field [H], we will not 
suppose a priori the order of the Gauss-Bonnet correction or the scalar field potential. 
Instead we will calculate leading-order gravity corrections for each of them, and obtain 
constraints on the relevant coupling constants (checking they fall within the validity of 
our perturbative expansion). Hence our analysis will apply for large couplings, which 
as we will see, are in accord with Gauss-Bonnet driven effective dark energy models. In 
this way we will show such models generally produce significant deviations from general 
relativity at local scales. We also include higher-order scalar field kinetic terms, although 
for the solutions we consider they turn out to be subdominant. 

In the next section we will present the theory in question and calculate the 
corrections to a post-Newtonian metric for a distributional point mass source. In 
section [3l we derive constraints from planetary motion, the Cassini probe, and a table- 
top experiment. For the Cassini constraint, we have to explicitly derive the predicted 
frequency shift for our theory, as it does not fall within the usual Parametrised Post- 
Newtonian (PPN) analysis. We discuss the implications of our results in section HI 

2. Quadratic Curvature Gravity 

We will consider a theory with the second-order gravitational Lagrangian 
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+ a [U(t^)CGB + 6(</>)G"^"V^(/.V,0 + e3(0)(V(/.)'VV + UmVct^Y] } , (1) 

which includes the Gauss-Bonnet term Cqb = - ^R^uR^"" + R^upaR^"'''' ■ Note for 
example that such a Lagrangian with given ,^'s arises naturally from higher dimensional 
compactification of a pure gravitational theory |15j. On its own, in four dimensions, the 
Gauss-Bonnet term does not contribute to the gravitational field equations. However 
we emphasise that when coupled to a scalar field (as above), it has a non-trivial effect. 

Throughout this paper we take the dimensionless couplings and their derivatives 
to be 0(1). There is then only one scale for the higher curvature part of the action, 
given by the parameter a, with dimensions of length squared. Similarly we assume that 
all derivatives of the potential V are of 0(V), which in our conventions has dimensions 
of inverse length squared. These two simplifying assumptions will hold for a wide range 
of theories, including those in which and V arise from a toroidal compactification of 
a higher dimensional space [15j. On the other hand it is perfectly conceivable that they 
do not apply for our universe, in which case the corresponding gravity theories will not 
be covered by the analysis in this article. 

Using the post-Newtonian limit, the metric for the solar system can be written [12] 

ds^ = -(1 - hoo){cdtf + {6ij + h,^)dx'dx^ + 0(6^/2) . (2) 

with /loo, hij = 0(e). The dimensionless parameter e is the typical gravitational strength, 
given by e = Gm/{rc^) where m is the typical mass scale and r the typical length scale 
(see below). For the solar system e is at most 10^^, while for cosmology, or close to the 
event horizon of a black hole, it is of order unity. The scale of planetary velocities v, 
is of order e^^^, and so the hoi components of the metric are 0(e^/2), as are dth^Q and 
dthij. In what follows, we will take (j) = (po + 0{e). For the linearised approximation we 
are using, we can adopt a post-Newtonian gauge in which the off-diagonal components 
of hij are zero. We can then write 

h,, = -2mij , hoo = -2$ , (3) 

and so 0^$ is the Newtonian potential. 

In this paper we will consider the leading-order corrections in e without assumptions 
on the magnitude of V and a. To leading order in e, the Einstein equations take the 
nice compact form. 



A$ = ^^Pm -V- 2<P($ + ^, 0) + 0(e^ ae'/r^ Ver^) (4) 



+ 0{e^aeyr^,Ver^) (5) 



^ , 47rG'o V 
= + Y - " 

where primes denote d/dcp, and V, etc. are evaluated at = (pQ. The matter energy 
density in the solar system is pm, and Go is its bare coupling strength (without quadratic 
gravity corrections). Other components of the energy-momentum tensor are higher order 
in e. The scalar field equation is 

A^ = V' -a [4e;P($, ^) + 6^^($ - ^, </>) + ^^'Dicp, (p)] + 0{e\ ae^r'', Ver^) . (6) 
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We have defined the operators 

AX = J2 X,u , V{X, Y)^J2 ^^^^^i - ^^^^ ■ 

i i,j 

with i,j = 1,2,3 where to leading order, the Gauss-Bonnet term is then Cqb — 
8r'($, For standard Einstein gravity {V — a — Q), the solution of the above 
equations is 

$ = * = (/) = (/)o, (8) 

where 

U„=*-^[t'.'^^. (9) 

J \x — x'\ 

We will now study solutions which are close to the post-Newtonian limit of general 
relativity, and take 

+ + = 00 + 50, (10) 

where 50, etc. are the leading-order a- and ^-dependent corrections. 

Note that the Laplacian carries a distribution and therefore we have to be careful 
with the implementation of the T> operator. We see that 50 is 0{V,a€^), and so, to 
leading order, we have 

A5(f)^V' -AaC[V{Um,Um). (11) 
Having calculated 50, we obtain 

A5$ = -V + 4<P([/^,50) (12) 
A5* = ^ + 2<P(C/„,50). (13) 
In the case of a spherical distributional source pm — mS^^\x), 

U^ = ^. (14) 

In accordance to our estimations for e the solar system Newtonian potentials are 
Um ^ 10~^, and the velocities satisfy v'^ < Um- For planets we have Um ^ 10~^ 
(with the maximum attained by Mercury). 
With the aid of the relation 

-n -m^ ^ '^nm ^ in+m+2) /^gX 

^ ' ' n+m+2 ^ ' 

the above expressions evaluate, at leading order, to 

6 c^r 



.+ ^-2C[^^^ (16) 



c^r 



c^r 



1 + 

1 + ^aV 



~6 7 c^r^ 



[17) 



rV 32{e,fa'iGomf 
+ 12""^ ^^^^ 
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We find that there are now non-standard corrections to the Newtonian potential which 
do not follow the usual parametrised expansion, in agreement with [19], but not [18] 
(which uses different assumptions on the form of the theory). First of all note that the 
Gauss-Bonnet coupling a couples to the running of the dark energy potential V , giving 
a 1/r contribution to the modified Newtonian potential (ITTll . We absorb this into the 
gravitational coupling. 



G — Gfi 



1 + ^-^aV 
3 



(19) 



The corresponding term in f|T8l) gives a constant contribution to the effective 7 PPN 
parameter. The r'^V terms in f[T7|) . f[T5]) are typical of a theory with a cosmological 
constant, whereas the final, 1/r'' terms are the leading pure Gauss-Bonnet correction, 
which is enhanced at small distances. If we take the usual expression for the PPN 
parameter 7 = \l'/$, we see that it is r dependent. In using the Cassini constraint on 7 
we must be careful to calculate the frequency shift from scratch. 

For the above derivation we have assumed 50 <C Um-, which implies V ^ Um/T"^ 
and a -C r"^ /Um- This will hold in the solar system if 

T . . „ -ifi 9 , I 10^^ m^ (everywhere) 

V < 10"^'' m-2 and « < < , ^,9 2 / 1 . ^ ^ 20 

I 10^ m^ (planets only) 

in geometrised units. Note that strictly speaking there is also a lower bound on our 
coupling constants, if the above analysis is to be valid. Indeed, if we were to find 
corrections of order ~ 10"^'', then it would imply that higher-order corrections from 
general relativity were just as important as the ones appearing in (fT7|) . (fT8l) . 

3. Constraints 

3.1. Planetary motion 

Deviations from the usual Newtonian potential will affect planetary motions, which 
provides a way of bounding them. This idea has been used to bound dark matter in the 
solar system ^20j, and also the value of the cosmological constant [21j. We will apply the 
same arguments to our theory. From the above gravitational potential (JT7j) . we obtain 
the Newtonian acceleration 



^d^ Gm 



Vr^ 64«)V2l Gmeff 



g(r) = -c^— = 1 - -- ^ 9 ^ _^™t ^21) 

where Vg = Gm/c^ is gravitational radius of the mass m. The above expression gives 
the effective mass rUes felt by a body at distance r. If the test body is a planet with 
semi-major axis a, we can use this formula at r a. Its mean motion n = a/ Gm/a^ 
will then be changed by 6n = {n/2){6mes/^)- By evaluating the statistical errors of 
the mean motions of the planets, 6n = —{3n/2)6a/a, we can derive a bound on 6mcs 
and hence our deviations from general relativity 

I 5m,s ^ _Va^ _ 64«)V| ^ 5a ^^2) 
3 m dvg 3a^ a 
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The values of a for the planets are determined using Kepler's third law, with a constant 
sun's mass niQ. Constraints on 5$ then follow from the errors 5a, in the measure of a. 
These can be found in [22], and are also listed in the appendix for convenience. Given 
their different r-dependence, the two corrections to 5mefr are unlikely to cancel. We will 
therefore bound them separately, giving constraints on a and V . 

The strongest bound on the combination ^'-^a comes from Mercury, with 

f^) <1.8xl0-^^ (23) 

Neglecting the cosmological constant term, and using a 5.8 x 10^ km and 1.5 km, 
we find 



< 



6rg 



^ 3.8 X 10^2 m^. (24) 

We see that this is within range of validity fl20|) for our perturbative treatment of gravity. 
In cosmology, the density fraction corresponding to the Gauss-Bonnet term is 



Qgb = 4^1"^^ • (25) 

If this is to play the role of dark energy in our universe, it needs to take, along with the 
contribution of the potential, a value around 0.7 at cosmological length scales (and for 
redshift 2; ~ 1). 

If we wish to accurately apply the bound on a (1241) to cosmological scales, details of 
the dynamical evolution of (p will be required. These will depend on the form of V and 
the ^i, and are expected to involve complex numerical analysis, all of which is beyond 
the scope of this work. Here we will instead assume that the cosmological value of 
is also 00) which, while crude, will allow us to estimate the significance of the above 
result. Given the hierarchy between cosmological and solar system scales it is natural 
to question this assumption but we will make it here, and discuss it in more detail in 
the concluding section. 

Making the further, and less controversial, assumption that dcp/dt ~ iJ, we obtain 
a very stringent constraint on VLgb'- 

I^gbI ^4|e;a|if2 <8.8x 10-=^°. (26) 

Hence we see that solar system constraints on Gauss-Bonnet fraction of the dark energy 
are potentially very significant, despite the fact that the Gauss-bonnet term is quadratic 
in curvature. 

Since we are assuming that all the are of the same order, the above bound 
also applies to the dark energy fractions arising from the final three terms in ([T]). 
Clearly there are effective dark energy models for which the analysis leading to the 
above bound ( l26l) does not apply. However any successful model will require a huge 
variation of ^1 between local and cosmological scales, or a very substantial violation of 
one of our other assumptions. 
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For comparison, we apply similar arguments to obtain a constraint on the potential. 
The strongest bound comes from the motion of Mars [21j, and is 



1^1 < 



1.2 X 10 



-40^-2 



(27) 



This suggests Vty = V/ (3-f^o) ~ '''■3 ^ 10^^, which is vastly weaker than the corresponding 
cosmological constraint {Vty ^ 1). Hence planetary orbits tell us little of significance 
about dark energy arising from a potential, in sharp contrast to the situation for Gauss- 
Bonnet dark energy. 



3.2. Cassini spacecraft 

The most stringent constraint on the PPN parameter 7 was obtained from the Cassini 
spacecraft in 2002 while on its way to Saturn. The signals between the spacecraft 
and the earth pass close to the sun, whose gravitational field produces a time delay. 
The smallest value of r on the light ray's path defines the impact parameter h. A 
small impact parameter maximises the light delay. During that year's superior solar 
conjunction the spacecraft was r^. = 8.43 AU = 1.26 x 10^^ m away from the sun, and 
the impact parameter dropped as low as 6min = 1.6i?0. A PPN analysis of the system 
produced the strong constraint 

57 = 7-1 = (2.1 ±2.3) X 10~\ (28) 

Given that our theory is not PPN we have to undertake the calculation from scratch. 

The above constraint comes from considering a round trip, in which the light travels 
from earth, grazes the sun's 'surface', reaches the spacecraft, and then returns by the 
same route. We take the path of the photon to be the straight line between the earth 
and the spacecraft, x = {x, b, 0) with x varying from —Xe to x®. For a round trip (there 
and back), the additional time delay for a light ray due to the gravitational field of the 
sun is then 



cAt = 

For the solution (fT7|) and 



dx = -2 



i^ + nr=V^dx. (29) 



181). this evaluates to 



cAt = Ar„ 



In 



452 



+ 



+ b (a® + re 



V 1024«; 
"6 ^ ¥~ 



'\2 3 

9 



(30) 



where we have assumed x® ~ a® ^ 6, and similarly for the spacecraft. 

Rather than directly measure At, the Cassini experiment actually found the 
frequency shift in the signal [23] 
_ dAt ^ dAt db 
~ ~dt db~di 
The results obtained were 

10-5 s db 



(31) 



b dt 



(2 + 57) 



(32) 
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If gravitation were to be described by the standard PPN formalism, then (57 would be 
the possible deviation of the PPN parameter 7 from the general relativity value of 1. 
From ( 130|) we obtain 

/ &V(ae + re) , 1536«)Vg Aai[V' \ Ar^db 
~ V 12r, + W 3 J cb df ^^^^ 



Requiring that the corrections are within the errors (l28l) of ( ]32l) . implies 

\^'a\<y^-<1.6xm^^m\ (34) 

This suggests the dark energy bound 

I^^gbI < 3.6 X 10-=^^ (35) 

although obtaining this bound from solar system data requires major assumptions about 
the cosmological behaviour of 0, as we will point out in section |H 

The data obtained by the spacecraft were actually for a range of impact parameters 
b, but we have just used the most conservative value b = femin = 1.6-R0. The above 
constraint is even stronger than ([23]), which was obtained for planetary motion. This 
is because the experiment involved smaller r, and so the possible Gauss-Bonnet effects 
were larger. 

Taking the above expression for y^^ ( l33l) at face value, we can also constrain the 
potential to be \V\ < 10~^^m^ and the cross-term |Q;,^^y| < 10~^. However these are 
of little interest as they are much weaker than the planetary motion constraints ([2] 
(!27|) . and also the former is far outside the range of validity (!20|) of our analysis. 



3.3. A table-top experiment 

Laboratory experiments can also be used to obtain bounds on deviations from Newton's 
law. For illustration we will consider the table-top experiment described in [23]. It 
consists of a 60 cm copper bar, suspended at its midpoint by a tungsten wire. Two 7.3 kg 
masses are placed on carts far (105 cm) from the bar, and another mass of m 43g is 
placed near (5 cm) to the side of bar. Moving the masses to the opposite sides of the bar 
changes in the torque felt by it. The experiment measures the torques iVio5 and —N^ 
produced respectively by the far and near masses. The masses and distances are chosen 
so that the two torques roughly cancel. The ratio R = Niq^/N^ is then determined, and 
compared with the theoretical value. The deviation from the Newtonian result is 

Sj, = - 1 = (1.2 ± 7) X 10"^ . (36) 

-'''Newton 

In fact, to help reduce errors, additional measurements were taken. To account for the 
gravitational field of the carts that the far masses sit on, the experiment was repeated 
with only the carts and a m' ~ 3 g near mass. The measured torque was then subtracted 
from the result for the loaded carts. 

The Gauss-Bonnet corrections to the Newton potential (fT7|l will alter the torques 
produced by all four masses, as well as the carts. Furthermore, since 5$ is non-linear 
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in mass, there will be further corrections coming from cross terms. The expressions 
derived in section [2] are just for the gravitational field of a single mass, and so will not 
fully describe the above table-top experiment. However, we find that the contribution 
from the mass m will dominate the other corrections, and so we can get a good estimate 
of the Gauss-Bonnet contribution to the ratio R by just considering m. 

The torque experienced by the copper bar, due to a point mass at X = (X, F, Z) is 



'bar Jhar ^ "'^ 



3 yX - xY 2 d^ 



(37) 

r=\X-x\ 



where pcu is the bar's density. A full list of parameters for the experiment is given 
in table I of [24]. The bar's dimensions are 60 cm x 1.5 cm x 0.65 cm. Working in 
coordinates with the origin at the centre of the bar, the mass m = 43.58 g is at 
X = (24.42, —4.77, —0.03) cm. Treating m as a point mass, Newtonian gravity implies 
a torque of ~ (8.2 cm^) Gmpcn is produced. The Gauss-Bonnet correction is 

QAGhrfiaC^ f , yX^xY ^ (Gm)3«)Vcu 
0X5 = pcn 2 / dx-^ -(0.025cm ) . (38) 

C Jbar \X — X\^ C 

To be consistent with the bound (l36l) . we require 6X5 /X^ to be within the range of Sr. 
This implies 

Kl < (18 cm^) ^ < 1.3 X 10^2 m^ (39) 

which is comparable to the planetary constraint (Ell). Extrapolating it to cosmological 
scales gives 

\nGB\<5.i X 10-=^°. (40) 

There are of course many more recent laboratory tests of gravity, and we expect 
that stronger constraints can be obtained from them. Table-top experiments frequently 
involve multiple gravitational sources, or gravitational fields which cannot reasonably be 
treated as point masses. A more detailed calculation than the one presented in section [2] 
will then be required. For example, the gravitational field inside a sphere or cylinder 
will not receive corrections of the form ( JTTll . and so any experiment involving a test 
mass moving in such a field requires a different analysis. 



4. Discussion 



We have shown that significant constraints on Gauss-Bonnet gravity can be derived from 
both solar system measurements and table-top laboratory experiments (note that further 
constraints arise when imposing theoretical constraints like absence of superluminal or 
ghost modes, see [25]). The fact that the corrections to Einstein gravity are second 
order in curvature suggests they will automatically be small. However this does not 
take into account the fact that the dimensionfull coupling of the Gauss-Bonnet term 
must be large if it is to have any hope of producing effective dark energy. Additional 
constraints will come from the perihelion precession of Mercury, although the linearised 
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analysis we have used is inadequate to determine this, and higher-order (in e) effects 
will need to be calculated. 

Performing an extrapolation of our results to cosmological scales suggests that the 
density fraction VLgb will be far too small to explain the accelerated expansion of our 
universe. This agrees with the conclusions of [19j. Hence if Gauss-Bonnet gravity is to 
be a viable dark energy candidate, one needs to find a loophole in the above arguments. 
This is not too difficult, and we will now turn to this question. 

In particular, we have assumed no spatial or temporal evolution of the field 
between cosmological and solar system scales, even though the supernova measurements 
correspond to a higher redshift and a far different typical distance scale. A varying cf) 
would of course imply that different values of ^j, and their derivatives, would be perceived 
by supernovas and the planets. It is interesting to note that the size of the bound we have 
found ( l26l) is of order the square of the ratio of the solar system and the cosmological 
horizon scales, s = (lAUifo)^ ~ 10~^°. Therefore one could reasonably argue that the 
small number appearing in (l26l) could in fact be due to the hierarchy scale, s, rather 
than a very stringent constraint on Qgb- This could perhaps be concretely realised with 
something similar to the chameleon effect [26] giving some constraint on the running of 
the quintessence theory. One other possibility is that the baryons (which make up the 
solar system) and dark matter (which is dominant at cosmological scales) have different 
couplings to (j) [22] • Again, this would alter the relation between local and cosmological 
constraints. 

Alternatively, it may be that our assumptions on the form of the theory should 
be changed. The scalar field could be coupled directly to the Einstein-Hilbert term, as 
in Brans-Dicke gravity. Additionally, the couplings and their derivatives could be of 
different orders. The same could be true of the potential. In particular, if were to have 
a significant mass, this would suppress the quadratic curvature effects, as they operate 
via the scalar field. This would be similar to the situation in scalar-tensor gravity with 
a potential, where the strong constraints on the theory can be avoided by giving the 
scalar a large mass (which, however, would inhibit acceleration). 

Finally, the behaviour of the scalar field could be radically different. We took it to 
be 0(e), like the metric perturbations. However since our constraints are on the metric, 
and not (f), this need not be true. Furthermore, since the theory is quadratic, there may 
well be alternative solutions of the field equations, and not just the one we studied. 

Hence to obtain a viable Gauss-Bonnet dark energy model, which is compatible 
with solar system constraints, at least one of the above assumptions must be broken. 
For many of the above ideas the higher-order scalar kinetic terms will play a significant 
role. This then opens up the possibility that the higher-gravity corrections will cancel 
each other, further weakening the constraints. We hope to address some of these issues 
in the near future. 
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Appendix 



For the benefit of readers without an astronomical background, we list relevant solar 
system parameters. The values for 5a come from table 4 of p2]. We take the Hubble 
constant to be Hq = 70kms~^ Mpc~^. 



Rq = 6.96 X 10** m 

r® = GniQ/c^ = 1477 m 

Ho/c = 7.566 X 10-27 

1 AU = am = 149597870691 m 



11 j-|^3 g 2 



G = 6.6742 X 10 
c = 299792458 ms-i 
tuq = 1.989 X 10^° kg 



kg 



-1 



name 


a (lO^m) 


6a (m) 


Mercury 


57.9 


0.105 


Venus 


108 


0.329 


Earth 


149 


0.146 


Mars 


228 


0.657 


Jupiter 


778 


639 


Saturn 


1433 


4.22 X 10^ 


Uranus 


2872 


3.85 X 10^ 


Neptune 


4495 


4.79 X 10^ 


Pluto 


5870 


3.46 X 10^ 
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